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Abstract 

In this paper, we develop a consistent extension of RPA for bosonic 
systems. In order to illustrate the method, we consider the case of the 
anharmonic oscillator. We compare our results with those obtained 
in mean-field and standard RPA approaches, with the exact ones and 
| show that they are very close to the exact ones. 
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1 Introduction 



X 

In many domains of physics perturbation theory fails to describe physical 
phenomena correctly. For example, bound states or collective excitations 
cannot be described perturbatively. In particular in the low-energy regime of 
strong interactions non-perturbative effects play an important role. Although 
much progress has been made in understanding the hadronic world, it is still 
indispensable to develop non-perturbative, symmetry conserving methods. 

In the last few years much effort has been made to apply non-perturbative 
techniques well established in many-body theory to quantum field theory 
(see, e.g., Refs. 0> @j B S ) ■ 111 ^ ms article we will concentrate on approaches 
related to the random phase approximation (RPA). In its standard form the 
RPA is known to respect symmetries. For example, in Ref. f|| it has been 
demonstrated within a linear cr-model that the RPA fluctuations generate 
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the Goldstone bosons related to the spontaneously broken global (chiral) 
symmetry. At mean-field level, on the contrary, these soft modes do not 
emerge. Thus the RPA seems to be a method very well suited to treat non- 
perturbative systems with spontaneously broken symmetries. Nevertheless 
it has some shortcomings: It is not of variational character and therefore 
does not necessarily fulfill the Rayleigh-Ritz criterion. In addition, it is not 
self-consistent, and therefore the intermediate states (corresponding to the 
mean-field quasiparticles) do not coincide with the final RPA particles. In 
particular, even if the RPA states show up as massless Goldstone modes, the 
intermediate states are massive. This has been pointed out by Aouissat et 
al. [0, U in the context of pionic excitations in an 0(iV)-model. 

It is obvious that this self-consistency problem can be circumvented in 
1 /iV-expansions [[3], [5], § . On the other hand, a fully self-consistent approxi- 
mation such as Self- Consistent RPA (SCRPA) M is desirable. This method 
has already been applied successfully to fermionic systems, e.g., to the Hub- 
bard model fH and the seniority model || . For a bosonic field theory a first 
attempt in this direction has been made in Ref. There some of the cor- 
relations are included, which are not present in the mean-field ground state 
used in standard RPA. This procedure, however, causes serious problems: 
In particular it does not preserve covariance, an essential ingredient of a 
relativistic quantum field theory. 

In order to establish a consistent method and to clearly state the problems 
related to the treatment of a bosonic system in an extended RPA formalism 
we will reduce the dimension of the system. Instead of a A$ 4 -theory as in 
Ref. || we will study a one-dimensional anharmonic oscillator. Of course this 
is a purely quantum mechanical, non-relativistic system and we can certainly 
not tackle the problem of covariance. Nevertheless, we can isolate some of 
the inherent contradictions related to bosonic systems. It has to be examined 
within a future work whether the method presented here is able to preserve 
covariance. 

Our toy model, the anharmonic oscillator has the advantage, that ex- 
act results are known and that we can directly compare the results ob- 
tained in different approximation schemes with those obtained by solving 
the Schrodinger equation. It has already been used as a test ground, e.g., for 
coupled cluster techniques |10[ or first attempts within SCRPA [p"T |. 



The article is organized as follows: We begin with presenting the RPA 
equations in Sec. |2|. In Sec. |3|, we review the Hartree-Fock (mean-field) ap- 
proximation together with the standard RPA. The renormalized RPA (r- 
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RPA) and its problems are discussed in Sec. [|. In Sec. we present our 
extended r-RPA approach. Sec. |6] is devoted to a discussion of the numerical 
results, in particular for the ground state energy, obtained within the different 
approximation schemes. Finally, in Sec. |7|, we will draw our conclusions. 



2 The RPA equations 

First of all, we have to state our starting point, the RPA equations. There 
are mainly two different ways to derive the RPA equations: the equation of 



motion method due to Rowe |12| and the Dyson-Schwinger-approach |13| . 
To derive the former, we suppose that we can generate an excited state \v) 
by the action of an operator Ql on the ground state |gs): |z/) = Qj,|gs), 
|gs) being the vacuum of Q u . Obviously Ql is a highly non-trivial many- 
body operator. Depending on the nature of |z/) (excited state with the same 
number of particles or excited state with a different number of particles), Ql 
can be seen as a superposition of (non-hermitian) one, two or more body 
operators. In general, Ql is expressed in the following form: 

Ql = J2KAl, (1) 

a 

where A^ a contains all the chosen excitation operators expressed in terms of 
<2j and a\ (annihilation and creation operators) and a represents all quantum 
numbers. Moreover, minimization of the energy E v = (u\B.\u) with respect 
to a variation 5Q U leads to a system of equations whose solution gives the 
amplitudes X v a and the excitation energies: 

(gs| [6Q V , [H, Ql]] |gs) = (E„ - E ) (gs| [6Q„, Ql] |gs) , (2) 

where E represents the energy of the ground state |gs). 

In the Dyson-Schwinger approach, which will be used in this paper, one 
starts from the definition of a time-ordered Green's function at zero temper- 
ature in equilibrium (the generalization to finite temperature is direct): 

G aP (t, = -<(gs|T (Aa(t) 4(f)) |gs) , (3) 

where T is the time ordering operator and A and A^ depend on time via 
A(t) = e tHt A{fd)e~ lHt , where H is the exact Hamiltonian. In equilibrium the 
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above two-time Green's function^ is a function of the time difference t — t' so 
that the Fourier transform G(E) depends only on one frequency. Our aim is 
to derive a Dyson equation for G(E). If we define 

G a p(t, t') = f ^e- iE WG a p(E), (4) 



the result is 



E G aP (E) = J\f a p + ^2 VC^M^ G ia (E) , (5) 



with 



Af af3 = (gs|K,Aj]|gs) , 
VC a(3 = (gs|[[A a ,H],4]|gs) . (6) 

It has to be noted that the operators entering the above Eq. (P) have to be 
taken at equal times since the truly dynamical contribution has been omitted 
(see, e.g., the derivation in Refs. @, @)- 

One can show that the two methods mentioned above are strictly equiv- 



alent provided that [14 



(gs|[H,[A,,4]]|gs> = 0. (7) 

Of course, if |gs) is the exact ground state, i.e., an eigenstates of the Hamil- 
tonian, this relation is automatically satisfied. Within an approximation this 
needs, however, not to be the case. Nevertheless, in most practical situations, 
including our treatment of the anharmonic oscillator, it can be shown that 
both methods are indeed equivalent. 

Obviously it is not possible to include the complete Hilbert space of ex- 
citation operators in A and it is necessary to restrict it. In principle the 
approximation can straightforwardly be refined by enlarging the space of ex- 
citation operators contained in A. But, in practice this often turns out to be 
a very difficult task because, on one hand, the dimension of the RPA matrix 
to be diagonalized increases and, on the other hand, there exist consistency 
problems for the theory itself (see the next section). 

1 Note that we deal with a two-time Green's function. That means that, if, e.g., A — 
a'<r, both particles are created at the same time. 
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Another crucial point in solving the above RPA equations, Eq. (||), is the 
way in which the expectation values are determined since the exact ground 
state is not known. It is however clear that it strongly influences the quality 
of the approximation. In standard RPA the exact ground state is usually 
replaced by a mean-field (Hartree-Fock) one. But we will show in the next 
section, using as an example the anharmonic oscillator, that it is possible to 
take into account an important part of the correlations present in the exact 
ground state in a fully consistent way. 



3 The model in standard RPA 

We start from the Hamiltonian of an anharmonic oscillator with quartic 
coupling: 

v 2 X 2 

K = P - + ^ + gX*. (8) 



with constants fx 2 < and g > 0. In the following sections, fi = y |/i 2 |. This 
Hamiltonian can be rewritten as 

r? Y 2 Y 2 

^ = P - + \ + ^-l)\ + gX' (9) 

and can be seen as a harmonic oscillator disturbed by a potential V(X) = 
(/i 2 — l)4r + gX A . We introduce the usual destruction and creation operators 
a and a' defined by a = (flX + ip) / y/2ji. The basis for standard RPA is the 
mean-field approximation which is usually formulated with the help of a 
Bogoliubov transformation. Thereby the starting point is the following trial 
wave function (see, e.g., Ref. [p!Qf| ): 

|HF) = e^ +1 ^ 2) \V) , (10) 

with variational parameters u and t. \V) is the vacuum of the initial destruc- 
tion operator a: a\V) = 0. The destruction operator of the Hartree-Fock 
vacuum 6|HF) = can be expressed with the help of u and t: 

b = (l-t 2 )- l l 2 (a-t(J -u) . (11) 

For later convenience we define to = — t)/{l + t) and s = ^2/jlu/(t — 1). 
In terms of the new creation and destruction operators we then obtain 

x = X + s = -L(6 + fot) ; 
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p = -^-(b-tf), (12) 

which makes the role of the variational parameters obvious: s simply de- 
scribes a translation in space and u a modified frequency of the oscillator. 
The Hamiltonian can be cast into the following form: 

Il = --(b-b^ + fi2 + 129S \ b + tff 
4 Auj 

^ + A9S \b + 6t) _ 4£L(6 + 6 t)s + J_ {b + 6 t)4 



2u V8^ 4w 2 



,2„2 



+^ 4 + ^f ■ (13) 

The aim of the Bogoliubov transformation is to select the best basis in the 
sense that the parameters u and s minimize the mean-field ground state 
energy, Eq F = (HF|H|HF). These two parameters are thus determined by: 

QfiHF dE HF 

q s = 0, ^ = 0. As can easily be checked by a direct calculation, the 
above equations are equivalent to the so-called "gap equations": 

8E HF 

-j^=0 * (HF|[H,6]|HF) = 0, 
8E HF 

--2— = <£> (HF|[H,66]|HF) = 0. (14) 

UUJ 

If |HF) was the exact ground state, the above gap equations would be satisfied 
automatically. However, the deeper origin of these equations (minimization 
of the ground state energy with respect to the Bogoliubov parameters ui and 
s) implies that they have always to be satisfied. 

Let us now discuss the RPA corrections to the ground state energy. 
In standard RPA the excitation operator A entering the RPA equations, 
Eqs. (H), is taken to be A 6 {b, tf, bb, b'tf}. We define the RPA ground state 
energy as follows: 

E^ FA = (RPA|H|RPA) . (15) 

This expression involves expectation values of the type (RPA|t4 5|RPA). 
These RPA corrected values can be calculated from the RPA Green's func- 
tions in the limit t — > t'\ 

r dE 

<RPA|A5|RPA> = % J i^G AB (E) . (16) 
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It is in this way that we introduce RPA fluctuations in the calculation of the 
ground state energy. 

In standard RPA the Green's functions on the right hand side of Eq. (|IED 
are not determined self-consistently, i.e., the expectation values entering the 
Green's functions themselves are determined using the mean-field ground 
state |HF). In our simple toy model we can obtain analytic expressions for 
the Green's functions. 

Note that from Eqs. (||) we do not directly obtain Green's functions like 
G x>x 2 or G PtP , for example, since the operators x, x 2 , and p are not contained 
in our set of excitation operators {b, b\ bb, tftf }. Hence, in order to apply 
Eq. (|16|) to expectation values like (x 3 ) or (p 2 ), the operators x and p are 
expressed in terms of b and (cf. Eq. (0)), and the resulting expectation 
values like (bbb), (b^bb), (bb), etc., are obtained from the Green's functions 
Gb,bb, G b i :bb , Gb t b, etc.0 The short-hand notation G x 3 or G p 2 has to be under- 
stood as the corresponding linear combination of Green's functions like Gb,bb, 
G b i bb , G b fi, etc., needed for the calculation of (x 3 ) or (p 2 ). 

To determine Ej} PA we need, in addition to (x 2 ), (p 2 ), and (x 3 ), also the 
expectation value (x A ). However, G x i contains Green's functions of the type 
G b i bbb , that is, A G {b, , bb, tftf} is not sufficient to calculate (x 4 ) []. We 
thus approximate this term by a factorization: 

(RPA|x 4 |RPA) ~ 3(RPA|a; 2 |RPA) 2 . (17) 

Finally, note that in the symmetric case (s = 0), the three-operator expecta- 
tion values vanish which implies that the Green's function are the Hartree- 
Fock ones: 

i i ^G X 7°(E) = (HF|x 2 |HF) (18) 
J Ztx 

Hence in this case the RPA energy is equal to the mean-field one. 

In the literature the term RPA is used for a similar approximation to the 
Green's functions as described above. But, it has to be mentioned that often 
the charging formula @, [15|] , 



E = E HF + / 1 ^(A|H mi (A)|A) , (19) 
Jo A 



2 Note that there is an ambiguity in the assignment of the time arguments to some of 
the Green's functions. This does, however, not influence the results (see Appendix ^ for 
a more detailed discussion). 

3 Since the excitation operators have to be non-hermitian, two-particle operators of the 
form b'b are not allowed. Thus we cannot use, e.g., G^tf, bb to determine G x i. 
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is applied to determine the ground state energy. Thereby H in4 (A) = XH int 
and the state |A) is defined as the ground state with respect to H(A) = 
H6 are + AHj ni where Hj„t is the interacting part of the hamiltonian. 

But, what we want to do is to compute directly (gs\H\gs) with \gs) = 
|HF), |RPA) or |0) (see next section) and compare the result obtained within 
the different approximations. Thus again, Eq PA is not the energy as under- 
stood in standard RPA but directly (RPA|H|RPA). 



4 The renormalized RPA approach 

What we want to do now is to go beyond the standard RPA approach and 
retain some of the correlations in the ground state. The starting point will 
be again the Bogoliubov transformation described in Sec. |3|. In contrast to 
standard RPA we will, however, not use the mean-field ground state |HF) 
but a more complicated one denoted |0). This new ground state, which will 
not be constructed explicitly in the following, contains some correlations. 
Since we want to use the "best" one, the parameters u and s will be again 
determined by the minimization of the ground state energy Eq = (0|H|0). 
This leads to generalized gap equations: 



dEo 

ds 
dE 



& (0|[H,6]|0) = 0, (20) 
<£> (0|[H,66]|0) = . (21) 



which are analogous to (0|[H, Qt]|0) = with the formalism introduced in 
Sec. |2|. We may note that the equivalence between these formulas is no longer 
valid if we consider three particle excitation operators (see Sec. [5] for further 
details). 

More precisely, independent of the approximation made to include corre- 
lations in the ground state, the above equations will guarantee that the basis 
is the energetically most favored one whithin this approximation. 

In the remaining part of this article we will assume that all expectation 
values of the type (0|66|0), (0|66^|0), (0|666^|0), etc. are real. This implies 
for example that (0|px|0) = —i/2 or (0|pa; 2 |0) = 0, independently of the 
specific ground state used. Thus, the key point is to determine the remaining 
expectation values in the a priori unknown ground state |0). In renormalized 



RPA [12], [IBL IT?], IT3 this is achieved with the use of Green's functions in the 
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limit t — > More precisely, we get self-consistency conditions for expectation 
values with two operators, i.e., (0|xa:|0) and (0\pp\0) from 

(0\AB\0) = i f^G AB (E) . (22) 

In contrast to Eq. ([16]), the expectation values entering the Green's functions 
on the right hand side of the above equation are determined self-consistently, 
i.e., not within the Hartree-Fock ground state. The remaining expectation 
values are approximated using a factorization, e.g., 

(0|pV|0) = (0|p 2 |0)(0|a: 2 |0) + 2(0\px\0) 2 . (23) 

This means that the self-consistent procedure is limited to expectation values 
with at most two operators and that the correlated part of the expectation 
values with more than two operators is neglected Q. Moreover, since (x) 
vanishes, it follows immediately that terms with three operators such as, 
e.g., (x 3 ) vanish. 

As for the standard RPA analytic expressions for the Green's functions 
can be given in terms of (x 2 ), (p 2 ) and (px) (see Appendix [A]). The latter 
expectation values have then, in principle, to be determined numerically via 
the self-consistency conditions, Eq. ( p2|) . With these values at hand, we 
are able, in principle, to calculate the ground state energy Eq. However, 
renormalized RPA is known to miss an important part of the correlation 



energy | 19fl . Besides, this method cannot be derived via a variational principle 



and therefore the Ritz criterion is no longer applicable. In fact, what happens 



in most cases, is that the ground state obtains overbinding |L9| . In our case it 
is even worse since one can show that a contradiction with the virial theorem 
and the Heisenberg uncertainty principle occurs. To that end let us look at 
the explicit equation for (p 2 ), which is obtained by integrating G p 2(E), given 
in Eq. (|35|), over E. The result is: 

(p 2 )= l J^G p2 (E) = e 2 (x 2 ) 



4(12sg) 2 (p 2 ) rdEl,- ^ 2w w> 
+t 1 J ^ ' J -^-( E ~ E !)(E 2 - E 2 2 ) , 

where e,Ei and E 2 are defined in Eq. (|36|) . If we now look at the explicit 
form (Eq. (|38D) of the gap equation for uj we conclude that the second term 
on the right hand side should vanish. However, this is only possible for {p 2 ) 
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= 0, implying a vanishing kinetic energy. This is obviously not a physically 
reasonable solution as it is in contradiction to the virial theorem. Even if we 
forget for the moment this contradiction, the self-consistent procedure can 
not give the right values for the occupation numbers and consequently for 
the energy. A possible way to circumvent that problem, that is to compute 
correctly the energy, is to use the "charging formula" as developed, e.g., in 
|| . It is applied in a similar way as in standard RPA (cf. Sec. |3|) to compute 
the correlation energy. The main difference is that we now have in principle 
a closed system of equations which have to be solved self-consistently. Fol- 
lowing [|J, one therefore replaces the gap equation concerning s, Eq. (pCf), 
by a minimization of the ground state energy obtained with the help of the 
charging formula (cf. Eq. ([H])). This has of course one disadvantage: One 
loses the symmetry of the double commutators, i.e. 

(lb, iK,bW}}) - ([&t 6 t ( [H ,6]]) = ([H,&t]> ^ o . (24) 

This symmetry property, however, is crucial for the stability of the sys- 
tem [Q. We will therefore proceed in another direction. 

5 The extended r-RPA approach 

Instead of using the charging formula, our strategy to include further correla- 
tions in a minimal way will be as follows: We will abandon the factorization 
approximation for the lowest non-vanishing expectation values of products 
of more than two operators. 

For the deformed case (s ^ 0) this means that we will keep the values of 
(x 3 ) and {p 2 x) (these are the only non- vanishing expectation values of three 
operators). Similar to (x 2 ) and (p 2 ), their actual values are determined via 
the self-consistency conditions, Eq. (|2"2"|). With these expectation values at 
hand the ground state energy is then obtained, as before, as the expectation 
value of the Hamilton operator (cf. Eq. flT3|)): 

E*° = ^ + ^((x 2 ) + s 2 ) + g(3(x 2 ) 2 + s 4 + 6s 2 (x 2 ) - As(x 3 )) . (25) 

First of all we have to note that there is no longer any obvious contradiction 
with the virial theorem. Moreover, as we will see in the next section, as far as 
we are not too close to the "phase transition" , our results for the ground state 
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energy are very close to the exact ones, without any further manipulation, 
like, e.g., the charging formula. This means (x 3 ) and {p 2 x) contain already 
an important part of the so far missing correlations. This does not necessarily 
imply that we also reproduce the exact values for, e.g., (x 2 ) or (p 2 ). But the 
correlations are combined in such a way that the energy is determined very 
precisely within this approximation. 

If we try to go further and include, e.g., also the correlated part of ex- 
pectation values with four operators, we face serious problems. To obtain 
the Green's functions necessary to calculate (x 2 ), (p 2 ), (px), (x 3 ) and (p 2 x) 
it suffices to restrict the space of excitation operators to A = {6, 6+ , 66, 6^6+}, 
i.e., two-particle excitations. As long as we stay in that space, the gap equa- 
tions together with the assumption of real expectation values for (66+), etc., 
guarantees the symmetry of the double commutators. As already discussed 
in Sec. [5], we need Green's functions of the type G b ^ bbb to obtain, e.g., (x A ). 
That is, three particle excitation operators have to included. In addition to 
the (practical) computational problems we encounter another difficulty: As 
mentioned in Sec. |3|, the Bogoliubov transformation contains only two pa- 
rameters and consequently, including three particles implies a priori an ap- 
proximation for the double commutators since there is no relation to enforce 
their symmetry. For instance, to ensure the following symmetry relation 

([bb, [H, 6+6+6+]]) - ([6+6+6+, [H, bb}]) = 6([H, 6+6+6]) , (26) 

we needed an equation which guarantees ([H, 6+6+6]) = 0. As mentioned 
above the Bogoliubov transformation does not provide us with an additional 
adjustable parameter. The assumption of a generalized gap equation of the 
form ([H, Q]) = as used in many RPA approaches also did not help. In this 
context boson expansion techniques |3], ||, mapping a pair of bosons onto one 
new boson, could be helpful. 

In the symmetric case s — 0, it is known that, by parity, (x 3 ) and (p 2 x) 
vanish. Following the same underlying philosophy as in the deformed case, 
we are now led to include (x 4 ), (p A ) and (p 2 x 2 ). This is, as stated above, a 
subtle task since the double commutators need no longer to be symmetric. 
Fortunately, this difficulty can be overcome: Although in principle there are 
inherent problems with the symmetry of the double commutators, one can 
explicitly check that all non-symmetric terms disappear for s = 0. This 
enables us to solve the Dyson equation without symmetrizing by hand (as it 
is often done in nuclear physics [pOfl). The corresponding Green's functions 
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are listed in Appendix [A|. The expectation values (x 2 ), (p 2 ), (x 4 ), (p A ) and 
(p 2 x 2 ) are then determined self-consistently using Eq. (|22|). In this case the 
ground state energy is given by: 

Er°= { -^ + ^(x 2 )+g(x*). (27) 

As can be seen in the next section, the numerical results for the ground state 
energy are in very good agreement with the exact ones. 

We have to emphasize that our approach is not fully self-consistent in 
the sense of SCRPA, although the expectation values of up to three or four 
operators, respectively, are calculated self-consistently. Let us mention one 
point which enlightens the difference of our method to the SCRPA approach 



used in Ref. |TT |. The authors use the following form for the RPA excitation 
operator: 

Q ] = -^(Afe^ 1 " -fxbb) . (28) 
V2 

This relation is easily inverted to give bb and frW in terms of and Q. Then 
the relation 

Q\SCRPA) = (29) 

can be applied to show that, e.g., (bb) vanishes. Our excitation operator 
contains in addition one and three particle operators, proportional to b, b' 
and bbbjtftftf , respectively. This renders the inversion less obvious and we 
therefore cannot directly make use of Eq. fl2"9"|). Thus, a priori, we have to 
calculate several expectation values which can be shown to vanish within 



the SCRPA formalism of Ref. |Tl[. The SCRPA results for the ground state 
energy in the symmetric case are also in very good agreement with the exact 
ones iTTl . 



6 Numerical results 

This section is devoted to a discussion of the numerical results. In Fig. H 
the energy of the ground state is displayed within the different approaches 
(exact result, mean-field, RPA, and our extended r-RPA result) as a function 
of n 2 /g. Of course, even the mean-field is in reasonable agreement with the 
exact results (10% at worst). The largest discrepancy can be observed in the 
region near the "phase transition" at /i 2 / ' g ~ —6.8, i.e., where s jumps from 
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some value s ^ to s = 0. This is also the case for the RPA and our r-RPA 
approach. But the figure clearly shows that the correlations incorporated in 
our approach are important: The agreement with the exact result is much 
better than in mean-field or RPA and discrepancies show up only in a very 
small region. Note that, the curve obtained within our approach is not 
continuous at the transition point. This is of course due to the difference 
in approximations made in the two domains: in one case, s ^ 0, one has 
incorporated the effect of (x 3 ) and (p 2 x), whereas in the other case, s = 0, 
we have included correlations in (x 4 ), (p 4 ) and (p 2 x 2 ). Near the transition, 
also higher order correlations become important, and we therefore miss the 
exact result near the transition point. 

Before we come to the conclusion, we would like to give one explicit 
example concerning the difference of our approach to the SCRPA: For fi 2 /g = 
— 1, we obtain (tfbbb) = —0.033 and (tftfbb) = 0.026, i.e., both expectation 
values are of the same order of magnitude. In the SCRPA approach used in 
Ref. |l ij the former expectation value vanishes identically. 

7 Conclusion 

We have considered an extended version of renormalized RPA that yields val- 
ues of high accuracy for the ground state energy of the anharmonic oscillator. 
Moreover, the self-consistent procedure developed here removes contradic- 
tions present in renormalized RPA, that is (x 2 ) and (p 2 ) can be determined 
without evidently violating the virial theorem or the Heisenberg uncertainty 
principle. 

This means that the correlations included, (x 3 ) and (p 2 x) for s ^ and 
(x 4 ), (p 4 ) and (p 2 x 2 ) for s = 0, represent the most important part of the 
correlations missed in renormalized RPA. Furthermore the systematics of 
the method allow us to think that it can be easily applied to quantum field 
theory and cure some of the problems which emerged in the work of Ref. Q . 
One problem of our approach, on which we certainly should improve, is the 
fact that we do not have a unified approach for the two cases s = and 
s ^ 0. 

Of course, it should be mentioned that there exist other powerful methods 
such as CCM (Coupled Cluster Method) [|l(J that describe very accurately 
the properties of the anharmonic oscillator. For problems concerning quan- 
tum chemistry or solid state physics, this technically and numerically very 



13 




o 1 

-15 -10 -5 

Figure 1: Exact solution for the ground state energy (solid line), mean field 
solution (dashed line), RPA solution (dotted line) and our extended r-RPA 
solution (dashed-dotted line) as a function of fi 2 /g with g = 0.5. Following 
[frnH , we added /i 4 /8 in order to have positive energies. In the region where 
s = the mean- field and the RPA solution coincide (see text). 



involved method and its extensions (see |T0[ for further details) are probably 
the most precise and well-suited. Actually, the approach developed in this 
paper probably can not compete with CCM in all these situations. Never- 
theless it has the great advantage to identify clearly the underlying physics 
(correlations) and to be very easy to use. 
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A Explicit expressions for the Green's func- 
tions 



We present here the analytical results for the various Green's functions 
needed throughout the calculation. Note that the Green's function listed 



One point is worth noting: We are dealing (in configuration space) with two- 
time Green's functions and the notation Ga,b{E) indicates that the time 
arguments for the corresponding Fourier transformed quantity Ga,b{^ — t') 
are assigned as A(t),B(t'). This means that there is an ambiguity in the 
assignment of time arguments to the Green's functions which contain more 
than two operators. In most cases this can be overcome by the fact that 
our excitation operators are limited to non-hermitian ones, i.e. an operator 
like tf(t)b(t) cannot exist. But, Gbb,bb(E) as well as Gb^bb(E) do exist and 
in general they do not lead to the same expression. We are, however, only 
interested in the limit t — *■ t' to determine the expectation values at equal 
times. In this limit all functions coincide. 

A.l Symmetric case 

In the symmetric case, i.e., s = 0, we obtain 



below do not directly emerge from the solution of the system of Eqs. ([5]) but 
we have to reexpress operators containing x and p with b and b*, cf. Eqs. ([12]). 
For instance 




GAE) 



(E 2 - £ 33 ) - Q % 



(E 2 — E 2 )(E 2 — El) 



G P 2{E) 



,2/ei2 v Q l2±gM; 




G :x *:(E) 



4a; 2 { w 3 (E 2 - E 01 E 02 ) 



(E 2 — Ef)(E 2 — El) 
1 / 2{X-2gu + ^{a 2 E + n)) 
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T 



G,AE) = -r ~ 



Aalu 6 (E 2 — E 2 )(E 2 — E 2 ) J ' 
uj 2 ( 2(A - 2gu + 3uj 3 (a 2 E + SI)) 



w^(E 2 - E 01 E 02 ) 

T \ 



G :p 2 x 2.(E) = 



4a 2 3 u 6 (E 2 - E 2 )(E 2 - El) 
a>2 E + E 02 



2uj 3 E 2 - E 01 E 2 ' 
where we have defined the following abbreviations: 



e = 

« 2 = 

a 3 = 

a = 

A = 



E 



CO 



J 33 — 

V = 

V = 

SI = 



T 



(31) 



/i 2 + 12^(a; 2 ) , 
(p 2 )+^ 2 (:r 2 ) 

l(--P 2 * 2 :)-3 + ^((: P 4 :) + 4(p 2 ) 

+Acu 3 (x 2 )+u 4 (:x 4 :)) , 
<p 2 )-u;V>, 

2#( - 3(p 2 )(l - 2uj(x 2 )) + uj(-2 
-6(: p 2 x 2 :) + 3co(5(x 2 ) - 2oj(x 2 ) 2 
+2cu(:x 4 :))))+ oj(uj 2 - e 2 )a , 

2 2 

1] Z — v A 

2 ' 

([666, [H,&W]]) , 
([666, [H,666]]) , 

^(\-4uJ 3 ((p 2 }+e 2 (x 2 )) + 22guJ 

+ 16w 2 #(-6(:r 2 ) + 3cj(:r 2 ) 2 - 2oj(: x 4 :))) , 

6a 3 (9A 2 (j + A (-27<j 2 + 4{a 3 E + rj + u)(E- u)) to 3 

+2au 5 (u (-4E(oi 3 E + rj) + 9cr 2 

+2(a 3 £ + 7] + z/V) + 2(a 3 £ + 77 - z/)e 2 ) ) , 
A — 2guj 



01,2 



a 2 



(X200 



(32) 
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The energies Ef and E\ are the solution of the following equation: 
E 2 

E 4 + (l8cr(A - au 3 ) - 4a 3 w 3 (£ 33 + e 2 )) 

+eV(36a 2 (r? - u) + 16a^S 33 )) = . (33) 

(a; 3 ) and (p 2 x) vanish because of parity. (: O :) denotes the expectation value 
of the normal ordered product of the operator O. If |HF) is taken as ground 
state those expectation values obviously vanish and (x 2 ) = l/(2u), (p 2 ) = 
u/2. From the gap equations, Eqs. fl2"U| , pT| ), we obtain the following condition 

uJ 2 ((p 2 } - e 2 (x 2 }) +3g- 12guj(x 2 }(l - lu(x 2 }) - Agu 2 (: x 4 :) = . (34) 

A. 2 Deformed case 

In the deformed case, i.e. s ^ 0, we obtain: 

GAE) 



(E 2 - E 2 )(E 2 - E 2 ) 



1 



F 2 Ji , 9g 2 «? 2 (/3- K ) 1 

16u 4 ti2 



GAE) = e 2 GAE) + ^ ^ " 



1 ' ~ x v ; a 2 (E 2 — Ef)(E 2 — El) ' 

68 g {a\E 2 + Aujk + 3E(k + a 2 uj)) 



G :X 3.(E) 



G :p 2 x: {E) = C12WE12 ^V" ' ( 6b > 



a 2 uj 2 (E 2 — E 2 )(E 2 — E. 
68g(a 2 E 2 - E(k + a 2 u)) 
a 2 {E 2 - E\)[E 2 - El 

where we have introduced in addition the following abbreviations 
e 2 = v 2 + 12g(s 2 + (x 2 )) 

P = A((p 2 ) + (x 2 )(e 2 + 12g(x 2 ))-18gs(:x 3 :)) 
8 = 2s(x 2 ) - (: x 3 :) 



2 ( (p 2 )(l + 4) " 12g 8 (: x 3 :) + (u 2 + e 2 )(x 2 ) 
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1,2 _ 



j) 2 + 8(12 S i) 



In this case we end up with two gap equations, one determining the value of 

s, 



Note that the expression for G x 2 (E) given in Eq. ([35|) can be interpreted as a 
free "propagator" with an additional "self-energy" depending on E. Similar 
interpretations are possible for all the other expressions. 



References 

[1] J.M. Hauser, W. Cassing, A. Peter, and M.H. Thoma, Z. Phys. A 353 

(1995) 301. 

[2] Z. Aouissat, G. Chanfray, P. Schuck, and J. Wambach, Nucl. Phys. A 603 

(1996) 458. 

[3] Z. Aouissat, P. Schuck, and J. Wambach, Nucl. Phys. A 618, 402 (1997). 

[4] H. Hansen, G. Chanfray, D. Davesne, P. Schuck, Eur. Phys. J. A14 (2002) 
397-411. 

[5] Z. Aouissat, Phys. Rev. C 62 (2000) 012201. 

[6] M. Oertel, M. Buballa, and J. Wambach, Phys. Atom. Nucl. 64 (2001) 
698. 

[7] G. Ropke, Z. Phys. B 99 (1995) 83; 

J. Dukelsky and P. Schuck, Phys. Lett. B 387 (1996) 233. 

[8] S. Schafer, P. Schuck, Phys. Rev. B 59 (1999) 1712. 




(37) 



and another one concerning u, 



(p 2 ) = e 2 (x 2 ) - 12gs(: x 3 :) . 



(38) 



18 



[9] A. Rabhi, P. Schuck, R. Bennaceur, G. Chanfray, and J. Dukelsky, 
preprint |nucl-th/ 0106064 

[10] R.F. Bishop and M.F. Flynn, Phys. Rev. A 38 (1988) 2211. 

[11] J. Dukelsky and P. Schuck, Nucl. Phys. A 512 (1990) 466; 
J. Dukelsky and P. Schuck, Mod. Phys. Lett. A 6 (1991) 2429. 

[12] D.J. Rowe, Rev. Mod. Phys. 40 (1968) 153 and Nuclear Collective Mo- 
tion, Methuen, London (1970). 

[13] P. Schuck and S. Ethofer, Nucl. Phys. A212 (1973) 269. 

[14] P. Schuck, J. Dukelsky, and G. Ropke, The Dyson equation approach to 
many body Green's functions, in Proceedings of the VI Hispalensis Inter- 
national Summer School, Oromana (Spain) June 1997, World Scientific, 
Singapore (1998). 

[15] A.L. Fetter, J.D. Walecka Quantum Theory of Many-Particle Systems, 
McGraw-Hill, New York (1971). 

[16] K. Hara, Progr. Theor. Phys. 32 (1964) 88. 

[17] F. Catara, N. Dinh Dang and M. Sambatoro, Nucl. Phys. A 579 (1994) 
1. 

[18] J. Toivanen and J. Suhonen, Phys. Rev. Lett 75 (1995) 410. 

[19] J. G. Hirsch, A. Mariano, J. Dukelsky, and P. Schuck Ann. Phys. 296 
(2002) 187. 

[20] M. Tohyama and P. Schuck, work in progress. 



19 



